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ABSTRACT 

The luminous Type Iln Supernova (SN) 2010jl shows strong evidence for the interac¬ 
tion of the SN ejecta with dense circumstellar material (CSM). We present observations 
of SN 2010jl for t ^ 900 d after its earliest detection, including a sequence of optical 
spectra ranging from t = 55 to 909 d. We also supplement our late time spectra and 
the photometric measurements in the literature with an additional epoch of new, late 
time BVRI photometry. Combining available photometric and spectroscopic data, 
we derive a semi-bolometric optical light curve and calculate a total radiated energy 
in the optical for SN 2010jl of ~ 3.5 x 10^^ erg. We also examine the evolution of 
the Ha emission line prohle in detail and hnd evidence for asymmetry in the prohle 
for t > 775 d that is not easily explained by any of the proposed scenarios for this 
fascinating event. Finally, we discuss the interpretations from the literature of the 
optical and near-infrared light curves, and propose that the most likely explanation 
of their evolution is the formation of new dust in the dense, pre-existing CSM wind 
after ^ 300 d. 
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1 INTRODUCTION 


Type Iln supernovae (SNe Iln) are a rare class of events, 
making up only 6 — 9 per cent of core collapse SN e (e.g., 
ISmartt et al.ll2009l : lLi et '^l201ll : ISmith et al]l2011^ b They 
are characterized by the presence of strong, composite- 
prohle emission lines of H and He. Generally, the emis¬ 
sion line prohles of SNe Iln may consist of narrow-width 
(NW, ~ 100 km s“^) a nd intermediate-width (IW, 
1000 km s“^) components (|Schlege]|ll99Ql : [FilippenkolllOOTl b 
The NW component is often explained as arising in dense, 
slowly expanding circumstellar material (CSM) photoion- 
ized by the initial flash of the supernova (SN) explosion 
fe.g.. IChevalier Fransson Il993: Salamanca et al.f 


Chugai et al.l I 2 OO 2 I : ISalamanca. Terlevich Tenorio-Tagld 


2003 ) • IW lines arise at early times as either Thomson 
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scatter broadening of the NW lines or emission from cool, 
accelerated, post-shock CSM gas (e.g. , IChugai &: Danzig^ 
I 994 I: I Chugai I 2 OOII : IChugai et al.l l2004 fpessart et al.ll2009l : 
Smith et al.l 201^ . The shock interaction of the ejecta 
with the dense CSM also provides a possible explana- 
tion for the extreme opt ical luminosity of many SNe Iln 
(jChugai Danzi^ll994l b 

SN 2010jl was discovered i n the galaxy UGC 5189A 
on 2010 November 3.52 (UT) bv lNewton &: PuckettI (|2010l b 
It was classihed as a SN Iln due to the presence of NW 
emission lines in a s pectrum tak e n on 2 010 November 5 
(jBenetti et al.l l2010l j . fsmith et al.l (|2011bl j identihed a pos¬ 
sible luminous blue progenitor star in archival Hubble Space 
Telescope (HST) WFPC2 data, but it is not possible to dis¬ 
tinguish between a single massive star and a star cluster 
from these observations. Observations of this region after 
the SN light fades could be used to differentiate between 
these possibilities. Regardless, each of the progenitor sce- 
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narios explored bv ISmith et, al.l feOllbh point to a progenitor 
with a mass of at least 30 M©. Spectropolarimetry results 
from SN 2010jl two weeks after its discovery indicate possi¬ 
ble asym metry in the explos ion geometry, or, more likely, in 
the CSM (|Patat et al.ll201lh . 

Pre-discovery images from the All-Sk y Automated Sur¬ 
vey (ASAS) North telescope in Hawaii (|Poimahskil l2002l : 
IPigulski et alJ l2009ll also show that SN 2010jl is intrinsi¬ 
cally luminou s, reaching a peak absolute /-band magnitude 
of ~ —20.5 (|Stoll et al.l 1201 ih . placing it near the class 
of < —21 mag ‘sup erluminous SNe’ defined bv iGal-YamI 
(|2012h . Additionally, IStoll et ^ (|201lh show that the host 
of SN 2010jl is low metallicity, supporting an emerging trend 
that optically luminous S Ne tend to go off in low-metallicity, 
low-luminosity host s fe.g.. lKozlowski et al.ll201^ : lNeill et al l 
l201ll:lLi et atll2nilh . The ASAS images indicate that the ex¬ 
plosion likely occurred in early October 2010. We will refer 
to days t since the earliest detection of S N 2010.il on 2010 
October 9 by ASAS (JD = 2,455,479.14; IStoll et al.llionh 
throughout this paper. 

[Andrews et al.l (|201lh construct a spectral energy distri¬ 
bution using optical, and near- and mid-infrared photometry 
of SN 2010jl at t ~ 90 d, and find evidence for a T ~ 750 K 
infrared (IR) excess. They interpret the IR excess as evi¬ 
dence for pre-existing du st in the unshocked CSM, heated 
by radiation from the SN. ISmith et al.l (|2012h observe a sys¬ 
tematic blueshift for the first seven months in the IW com¬ 
ponent of Ha that arises in the post-shock CSM, indicating 
the formation of new dust in the post-shock cooling zone 
of the SN. Their spectra also indicate that the blueshift is 
stronger at shorter wavelengths, which one would expect if 
it is the result of dust formation; however, this explanation 
co nflicts with the low d ust emission temperature observed 
bv [Andrews et ahl (|201lh . An alternative explanation is that 
the high optical depth continuum continues to block light 
from the receding side of the post-s hock shell, resultin g in 
line profiles that appear blueshifted (jSmith et al.[[2Ql3 l. 

An extensive data set consisting of UBVRI light curves 
and a seq uence of opti c al spe ctra for over 500 days was pre¬ 
sented in [Zhang et al.[ (|2012[ ) . Both the slow decline of the 
light curves and the remarkable strength of the Ha emission 
with respect to the continuum for over 500 days provide 
strong evidence for successive interactions of the SN ejecta 
with a dense, hydrogen-rich CSM. They also observe a sys¬ 
tematic blueshift of the line profiles of Ha, but disfavor the 
interpretation of new dust formation due to a decrease in the 
blueshift of the line profiles for t > 500 d, and the lack of a 
corresponding decline in the optical light curves one would 
expect from extinction by new dust. 

[Maeda et al.[ (|2013[ j claim to find strong evidence for 
dust formation in SN 2010jl at t ~ 550 d. They present near- 
IR (NIR) spectroscopy and JHK photometry, in addition 
to optical spectroscopy and BVRI photometry at similar 
epochs. They detect a NIR thermal component with tem¬ 
peratures of T ~ 1000 — 2000 K which was not present at 
earlier epochs. They also report an increase in the blueshift 
of Ha to > 700 km for t > 400 d, and that the degree of 
the blueshift is smaller for lines at longer wavelengths. Fi¬ 
nally, they point out that the decline in optical luminosity 
accelerates after ~ 1 yr. They conclude that the cumulative 
evidence provides a strong case for new dust formation in 


a dense cooling shell formed by the interaction of the SN 
ejecta and the CSM. 

A more complete data set including uBVRiJHK pho¬ 
tometry, ultraviolet (UV) spectra from HST, extensive op- 
tic al spectroscopy , and spectra in the NIR was presented 
in [Fransson et al J ([2014[j. in which t he au thors debate the 
conclusions drawn bv iMaeda et al ] (l2013h . They derive a 
‘pseudo-bolometric’ (optical through NIR) light curve of 
SN 2010jl, and show that the NIR component, which they in¬ 
terpret as an echo from pre-existing CSM dust heated by the 
SN, becomes an increasingly important contribution to the 
bolometric luminosity with time. They present a sequence 
of Ha profiles from t = 31tot = 1128d, and claim that 
the broad component shows smooth, peaked, and symmet¬ 
ric profiles characteristic of electron scattering. The blueshift 
present in the lines for t > 50 d is interpreted as a macro¬ 
scopic bulk velocity due to radiative acceleration of the pre¬ 
shock scattering medium by the energetic radiation field of 
the SN. From the narrow CSM lines present in their spectra, 
some of which show a P Cygni absorption component, they 
infer a CSM velocity of ~ 100 km s“^ consistent with a Lu¬ 
minous Blue Variable (LBV) star as a possible progenitor. 
They find the mass loss rate to account for the bolometric 
light curve with CSM interaction to be > 0.1 M© yr“^, and 
a total mass lost of > 3 M©. 

[Borish et al.[ (|2015[j present a sequence of 1 — 2.4 /xm 
NIR spectra from t = 36 to 565 d. They observe a progressive 
blueshift to ~ 700 km s“^ of the broad hydrogen lines with 
time, similar to the evolution seen in the optical hydrogen 
lines. They also observe a broad, asymmetric component in 
the He i AA10830 and 20587 lines, and a blue shoulder feature 
at ~ 5000 km s“^, possibly associated with a high velocity, 
He-ric h ejecta flow. Similar to the results of [Fransson et aH 
they observe a T ~ 1900 K blackbody feature in the 
continuum emission at t = 403 d, likely associated with dust 
emission. 

iGall et all (l2014l^ observed SN 2010jl with the VLT/X- 


shooter spectrograph in a sequence of optical-NIR spectra 
covering the range 0.3 — 2.5 /xm for t < 250 d, and again 
on t = 876 d. They find evidence for wavelength-dependent 
extinction in their analysis of the IW emission lines, and in¬ 
terpret this as extinction by newly formed dust. They derive 
dust extinction curves for SN 2010jl and find evidence for 
large dust grains in excess of 1 /xm in size. They claim that 
the dust production site transitions from the CSM to the 
ejecta for t > 500 d based on the accelerated growth in dust 
mass at this time. 

SN 2010jl w as detected by the S wift on-board X-ray 
Telescope (XRT, [Burrows et al.[ [2005[ i on 2010 November 
5.0 — 5.8 with an inferred X-ray luminosity of (3.6 ± 0.5) x 


10^° erg s“^ (jlmmler. Milne ^ Foo\^\2QlQ\ Chan dra ob¬ 
servations were triggered by Chandra et al. (|2012a[ i to ob¬ 
serve SN 2010jl on t = 59 and 373 d, which demonstrated a 
high temperature of > 10 keV for both epochs, and a high 

absorption column density of ~ 10^^ cm“^ at the first epoch 
that decreases by a factor of 3 at the second epoch, likely 
asso ciated with absorp tion by the CSM. 

[Ofek et m] (|2014[ i also present X-ray observations of 
SN 2010jl from NuSTAR and XMM-Newton, in addition to 
multi-epoch Swift-XRT observation s, and public Chandr a 
observations (Pis Pooley, Chandra: [Chandra et al.[[2012b[ i. 
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They also present the PTF R- and ^-band, and Swift-VYOT 
light curves. They apply a more general version of the model 
described bv ISvirski et al.l (|2012l l for an SN shock interacting 
with an optically thick CSM, assuming spherical symmetry, 
and hnd a CSM mass in excess of 10 Mq. Furthermore, they 
measure the temperature of the X-ray emission and infer a 
shock velocity of ~ 3000 km s“^. 

In this paper, we present an analysis of extensive op¬ 
tical observations of SN 2010jl over a period of ~ 2.5 yr. 
Our data set adds significant coverage of the spectroscopic 
time series for SN 2010jl and we evaluate the the current 
physical models proposed in the literature to explain these 
observations. 

The observations, primarily a sequence of optical spec¬ 
tra but including a supplementary epoch of late time BVRI 
photometry, and data reduction procedures are presented in 
Section O In Section O we describe our analysis of these 
data, including the generation of a semi-bolometric, opti¬ 
cal light curve in Section 13.2[ and a detailed study of the 
evolution of emission line prohles, particularly Ha, in Sec¬ 
tion 13.31 Finally, we discuss our physical interpretation of 
these results and comment on previous interpretations in 
Section [H We summarize our conclusions in Section O 


2 OBSERVATIONS 

In this section, we present extensive optical spectroscopy of 
SN 20I0jl and an additional epoch of late time photometry, 
and describe data reduction procedures. 

2.1 Optical spectroscopy 

We obtained a total of twelve optical spectra of SN 20I0jl us¬ 
ing a number of ground based instruments: the Dual Imaging 
Spectrograph (DIS) on the 3.5-m ARC telescope at Apache 
Point Observatory (APO), the Wide Field Reimaging CCD 
Camera (WFCCD) on the du Pont Telescope at Las Cam- 
panas Observator y (LCO), the Mu lti-Object Double Spec¬ 
trograph fMODS JPogge et al.ll20IQl l on the Large Binocular 
Telescope (LBT) at Mt. Graham International Observatory 
(MGIO), and the Ohio State Multi-Object Spectrograph 
fOSMOS. IStoll et~aLll20I0l : lMartini et al.ll20Ilh on the 2.4-m 
Hiltner Telescope at MDM Observatory. These observations 
began on t = 55 d and continue until t = 909 d. These 
spectroscopic observations are summarized in Table [U Also 
listed in Tabled] ar e the two spectros copic observations orig¬ 
inally presented in IStoll et ah] (|20III L which are used in the 
analysis discussed below in Sections 13.31 and [3.3.II 

The Hiltner/OSMOS spectrum, the APO/DIS spectra, 
and the du Pont/WFCCD spectra were reduced using stan¬ 
dard tasks in iRAij^. The LBT/MODS spectrum was reduced 
using routines for MODS dual-channel long-slit grating spec¬ 
tra following the ‘MODS Basic CCD Reduction with mod- 
sCCDRed^ manuaQ. Cosmic rays wer e removed from th e 
images for all spectra using LACOSMIC (jvan Dokkumll200~lh . 

^ IRAF is distributed by the National Optical Astronomy Observa¬ 
tory, which is operated by the Association of Universities for Re¬ 
search in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 

^ modsCCDRed manual is available here: 


Table 1. Spectroscopic observations of SN 20I0jl 


UT Date 

Phase 

(d) 

Range 

(A) 

Instrument 

2010 Nov. 6 

28 

3920-6800 

OSMOS 2.4 m Hiltner 

2010 Nov. 12 

34 

3110-5880 

OSMOS 2.4 m Hiltner 

2010 Dec. 3 

55 

3460-9500 

DIS 3.5 m ARC 

2010 Dec. 29 

81 

3660-9100 

WFCCD 2.5 m du Pont 

2011 Mar. 12 

154 

3460-9400 

DIS 3.5 m ARC 

2011 Nov. 17 

404 

3170-9890 

MODS LBT 

2012 Jan. 02 

450 

3460-9500 

DIS 3.5 m ARC 

2012 Jan. 20 

468 

3460-9700 

DIS 3.5 m ARC 

2012 Feb. 22 

501 

3460-9700 

DIS 3.5 m ARC 

2012 Mar. 4 

512 

3830-7020 

DIS 3.5 m ARC 

2012 Nov. 22 

775 

3920-6800 

OSMOS 2.4 m Hiltner 

2013 Jan. 29 

824 

3660-9100 

WFCCD 2.5 m du Pont 

2013 Mar. 18 

891 

3460-9500 

DIS 3.5 m ARC 

2013 Apr. 6 

909 

3660-9100 

WFCCD 2.5 m du Pont 


The phase is the number of days since earliest detection on JD 
= 2,455,47 9 .14. T he first two spectra were originally presented in 

IStoll et alJ ll201lh . 


Wavelength and flux calibrations were performed with com¬ 
parison lamps and standard star spectra, respectively. The 
small wavelength zero-point shift between the comparison 
lamp wavelengths and the observed wavelengths in the spec¬ 
tra was also corrected. 

The spectra were Doppler corrected for the redshift of 
the host galaxy [z = 0.0107, NEeB) and for Galactic red¬ 
dening using E{B — U) = 0.024 from ISchlaflv fc Finkbein^ 
(|20Ilh . Following ISmith et al.l (l20IIblL no co r rectio n was 
made for host galaxy extinction. iFransson et al ]<|^ made 
an estimate of the reddening for both the host galaxy and the 
Milky Way from the damping wings of the Lya absorption in 
their spectra, and they adopted a value of E{B — V) = 0.058 
for the total reddening. Our spectra are shown in Fig. [T] 
showing the log of the flux plus an arbitrary constant for 
clarit y. We also show the two early time OSMOS spectra 
from istoll et al.l (|20Ilh . 

2.2 Optical photometry 

In order to better interpret the spectra, we performed flux 
calibr ations (see Sectio n 13.21) using publi s hed p hotometry 
from IStoll et m] (|20Ilh and IZhang et al.l (|20I^ ). and ob¬ 
tained an additional epoch of late time BVRI imaging of 
SN 20I0jl using OSMOS on the 2.4-m Hiltner Telescope 
at MDM Observatory on t = 776 d. The images were re¬ 
duced using standard tasks in the IRAF ccdproc package. 
Flat-fielding was performed using twilight sky flats for each 
hlter taken from t = 773 to t = 776 d. The photometry 
of the SN was obtained using the digital p hotometry re¬ 
ducti on program DAOPHOT ll with ALLSTAR (|Stetsonlll987l . 
l2000li . Due to issues with over subtracting the source, we 

http://WWW.astronomy.ohio-state.edu/ 
MODS/Manuals/MODSCCDRed.pdf 

^ The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet propulsion Laboratory, California Institute of Tech¬ 
nology, under contract with the National Aeronautics and Space 
Administration. 
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Rest Wavelength (A) 

Figure 1. Rest wavelength spectra of SN 2010jl obtained 
from November 2010 to April 2013. The spectra are cor- 
rected for Galactic extincti on with E(B — V) = 0.024 mag 
(jSchlaflv fc Finkbeinerl [SoTlh . Each spectrum is shifted by the 
indicated constant for clarity. Days since first detection are also 
indicated along the right side. A lso shown at the t op are the two 
early time OSMOS spectra from IStoll et al.l (l201lh . 


Table 2. Photometry of SN 2010jl from JD = 2, 456, 255 


Filter 

m 

(mag) 

(7 

(mag) 

M 

(mag) 

fj 

(mag) 

(Lq) 

B 

19.31 

0.11 

-14.22 

0.11 

4.3 X 10'^ 

V 

18.94 

0.05 

-14.57 

0.05 

4.1 X 10'^ 

R 

17.13 

0.05 

-16.36 

0.05 

1.6 X 10® 

I 

17.87 

0.05 

-15.60 

0.05 

4.9 X 10^ 


tested various annuli in DAOPHOT to estimate the back¬ 
ground contribution and to minimize contamination from 
the host galaxy. Instrumental magnitudes were converted 
to ap parent magnitudes usin g nearby field stars with SDSS 
DR7 (jAbazaiian et al.l l2009l i ugriz ph otometry, converte d 
to standard BVRI magnitudes using llvezic et al.l (|2007f) . 
Our photometric results are summarized in Table [21 where 
apparent magnitudes are converted to absolute magnitudes 
by adopting a distance modulus to the host galaxy of 
33.43 mag (NED) and using the Galact ic extinction correc¬ 
tions from ISchlafiv fc Finkbein^ (|201ll ). BVRI magnitudes 
are also converted to band luminosities, vL^ILq^ using 
the zero poin t s of t he Johnson-Cousins-Glass system from 
iBessell et al.l (|l998l ). Our photometric results are show n 
in Fig. [21 al o ng w it h light curves fro m IStoll et ^ (1201 ll), 
IZhang ^ al.l (|2012l ) , lOfek et ^ (|2014l ) , and Fransson et al l 



Figure 2. BVRI light curves of SN 2010jl. The stars are the new 
data presented in this p aper. Circles are early time measurements 
from ^oll et all (l20ll1v The diamonds c ome fro ny[Zha^ et^l 
(|2Q12 b B- and E-band measurements from iFransson e'^dT(|^14l ) 
are shown as tri angles. The ‘X’ sym bols are the PTE R-band mea¬ 
surements from lOfek et (|2014l i , co nverted from AB to Veg a 
magnitudes using the conversions from iBlanton &: RoweisI (|2007l b 
Error bars are shown, but can be smaller than the points. 


3 ANALYSIS 

In this section we present the results of our data analy¬ 
sis. In Section nm we examine the new photometric data 
points presented in this paper in the context of published 
light curves from the literature. In Section [321 we describe 
the procedure used to derive a semi-bolometric optical light 
curve for SN 2010jl, and compare its evolution to that of Ha 
and the NIR luminosity. In Section [T3l we describe the evo¬ 
lution of the optical spectra, in particular the Ha emission 
profile (Section 13.3.II) . 


3.1 BVRI light curves and color evolution 


The photometric measurements from 2012 November 23 pre¬ 
sented in this paper are added to light curves from the lit¬ 
erature as shown in Fig. [21 The light curves also include the 
pre-di scovery V- and /-band measu rements fronilStoll et al.l 
boilh . the BVRI light curves fro ml^ang et aP (l2012l ). the 
PTF R-band measurements from lOfek et al.l (| 20 111 ), con- 
ve rted from AB to Vega ma gnitudes using the conversions 
by [Blanton Ro weisI (l2 007l) . and the B- and E-band light 
curves from iFransson et al.l (|2014l) . We note a few discrepan¬ 
cies among the measure ments, particularly at late times. Be¬ 
ginning at t ~ 80 d, the [Fransson et HI (l2014l ) B and E light 
curves begin to diverge fro m the corresponding Zhang et al.l 
(|2012l ) light curves, with the lFransson et al ](Eil measure¬ 
ments ~ 0.1 mag fainter in B and ~ 0.2 mag fainter in E. Af¬ 
ter the break in the light curve at t ~ 300 d, t hese discrepan¬ 
cies gr ow to ~ 1.0 mag, and we note that the lFransson et aH 
(|2014l ) light curves show a significantly steeper decline in 
both bands. Our measurements at t = 776 d, shown as stars 
in Fig. [21 are ~ 0.2 mag fainter in R than the PTF mea- 
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Figure 3. Color evolution of SN 2010jl. The green stars were ob¬ 
tained using the BVRI p hotometry presented in this paper. The 
blue diamonds co me from IZhan g et al.l (I 2012 I) . and the cyan tri¬ 
angles come from iFransson et al.l (|20l3 b All measurements were 
corr ected for Galactic reddening using the extinction corrections 
from ISchlaflv &: Finkbeineil (1201 ll ) before computing colors. Error 
bars are shown, but can be smaller than the points. 
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sur ements and ^ 0 . 2 — 0 .3 mag brighter in B and V than 
the IFransson et ahl (|2014l l measurements. These discrepan¬ 
cies may be due to difficulties in estimating the background 
flux as the SN fades and the contribution of the host galaxy 
becomes more significant. 

Fig. E] shows color measurements from our BVRI 

photometry as _ compa red to the co lor _ ^olution from 

IZhang et al. I (l2012h and IFransson et 'Zl (l2014h . The B - V 
color remains mostly constant at ~ 0.3 mag, consistent with 
our measurement of B — V = 0.3 zb 0.1 mag at t = 776 d. 
In B — V , we again no te a ^0.1 mag discrepanc y between 
the IZhang et al.1 (|2012l i and IFransson et al.1 (|2014h measure¬ 
ments. The B — V color al so appears to fa l l off t oward the 
blue for t > 800 d in the IFransson et al. I Jioli) measure¬ 
ments, which may be associated with the bluer continua 
present in the last few spectra, possibly from host galaxy 
contamination. In V — R, the color becomes redder with 
time, increasing from V — R ^ 0.4 mag at early times to 
our measurement of V — R = 1.8 zb 0.1 mag at t = 776 d. 
This is consistent with the continued strength of the Ha line 
(contained in the i?-band) relative to the continuum at late 
times (see Sections 13.31 and 13. 3. ID . We also note that V — I 
becomes redder, increasing from V — I ~ 0.6 mag at early 
times to our late time measurement ofV — I= l.OzbO.l mag. 


3.2 Semi-bolometric optical light curve 

We generated a semi-bolometric, optical light curve (3460 
- 9700 A) using a combination of the spectra and avail ¬ 
able photometric data from this paper, IStoll et ^ (|201lh 
and IZhan g et al. ( 2012li. We also used the 12 optical spectra 
from IZhang et ah (|2012l i. taken from the online Wei zmann 
Interactive Supernova data REPository (WISeREP, lYaronI 


I 2 OI 2 I 1 FI. The shape of the optical spectral energy distribu¬ 
tion (SED) was estimated as a function of time by interpo¬ 
lating the observed spectra, put on an absolute flux scale 
using synthetic photometry. Using this model of the SED 
with time, synthetic spectra were produced for days that 
had available photometry, and integrated to give the op¬ 
tical luminosity for each day. The resulting light curve is 
tabulated in Table [3] and shown in Eig. |4l Statistical un¬ 
certainties were estimated for each point using a bootstrap 
technique. The entire analysis was run 10,000 times, allow¬ 
ing the input photometry to vary with Gaussian distributed 
noise defined by the reported uncertainties. Trapezoidal in¬ 
tegration of this light curve in time from t = 0 to 776 d 
gives an estimated total radiated energy in the optical of 
(3.50 zb 0.03) X 10^° erg. We caution that the uncertainty 
reported here is only an estimate of the statistical uncer¬ 
tainty from the bootstrap procedure, and does not include 
systematic error such as host galaxy contamination and the 
gap in the data when the SN was too close to the Sun to be 
observed between t ^ 200 and 40(3 d. Th is number is similar 
to that derived bv lEransson et aP (|2014l i of ~ 3.3 x 10^° erg 
in the 3600 — 9000 A range from t = 0 to 920 d. In addition 
to the difference in wavelength and tim e coverage, w e point 
out again the discrepancies between the IZhang et al.1 (|2012l i 
photometry used i n this analysis and the photometric mea¬ 
surements used by lEransson et aP (|2014h (see Section IQ) , 
which may be responsible for the difference in the derived 
optical energy output. 

Overplotted as black dashed lines in Eig. 0] are power 
law fits to the derived optical light curve. The points from 
t ~ 20 to ~ 370 d were best lit with a power law decay, 
given by Lopt(t) = 1.02 x 10^^(t/100 erg s“^, where 

Lopt is the optical luminosity and t is the time since the 
earliest detection of this event. The remaining late time 
points through t — 776 d were fit with an Lopt(t) = 
6.01 X 10^^(t/370 erg s“^ power law. The break in 

the light curve, defined as the intersection of the two power 
law fits, occurs at thr ~ 370 d after first detection. The lack 
of data near this epoch makes it difficult to accurately con¬ 
strain the timing of the break. Possible interpretations of the 
decline in the optical light curve at late times are discussed 
in Section 01 

Eor comparison in Eig.|4l we also show the Ha luminos¬ 
ity as a function of time (also see Eig. [5] and Section [ffi3l for 
further discussion) and the corresponding power law fits. 
Luminosities corresponding to open triangles are only ap¬ 
proximate due to insufficient wavelength coverage in the 
spectra for proper flux calibration with synthetic photom¬ 
etry and are not included in the power law fits described 
below. In contrast to the integrated optical luminosity, the 
early time Ha luminosity follows an increasing power law 
given by Luait) = 7.94 x lO'^^Ci/lOO erg s“^ un- 

til at least t = 154 d. At late times, the Ha light curve 
is well fit by steep power law decay given by Luait) = 
1.43 X 10^^(t/350 d)“^'^® erg s“^. The slope of the Ha light 
curve remains slightly shallower than that of the optical light 
curve at late times. The break in the Ha light curve, defined 


^ WISeREP spectra are available here: 
http://wiserep.weizmann.ac.il/ 


© 0000 RAS, MNRAS 000 , 000-000 


































































6 J. E. Jencson et al. 


Table 3. Optical luminosity 


JD Phase L L 

(d) (ergs-i) (Lq) 


55479.14 

0.00 

1.6(2) X 

10"^® 

4.2(4) X 10^ 

55484.13 

4.99 

i.ej) X 

10« 

4.1(3) X 10^ 

55488.12 

8.98 

I.tJ) X 

10"^® 

4.4(3) X 10^ 

55493.14 

14.00 

1.75) X 

10"^® 

4.4(4) X 10^ 

55494.12 

14.98 

1.7(2) X 

10^3 

4.4(4) X 10^ 

55504.10 

24.96 

i.ej) X 

10^3 

4.1(3) X 10^ 

55509.10 

29.96 

l.sj) X 

10« 

3.9(4) X 10^ 

55509.13 

29.99 

l.lj) X 

10"^® 

3.7(3) X 10^ 

55515.75 

36.61 

1.565) X 

10« 

4.01(9) X 10^ 

55516.75 

37.61 

1.585) X 

10"^® 

4.0(1) X 10^ 


The phase is the number of days since earliest detection on JD = 
2,455,479.14. Table [3] is published in its entirety in the electronic 
edition of this journal. A portion is shown here for guidance re¬ 
garding its form and content. 



Figure 4. The semi-bolometric optical light (‘X’ symbols), Ha 
(triangles), and AT-band luminosities (squares) of SN 2010jl. The 
error bars, representing only statistical uncertainties for the de¬ 
rived optical points, are shown but can be smaller than the points. 
Power law fits to the optical light curve and Ha are shown as 
black dashed and blue dotted lines. The Ha luminosities associ¬ 
ated with open triangles are approximate and are not included in 
the fits. 


by the break in the power law fits, occurs at t ~ 350 d, near 
the time of the break in the full optical light curve. 

Finally, w e als o show the AT-band data from 
iFransson et al ] (|2014l ) in Fig. [H converted from 2MASS 
system magnitudes to lyLjy band luminosities using a Ks 
hux zero point of — 666.8 Jy at u = 1.390 x 

10^^ Hz (ICohem Wheaton fc MegeathI 1200 j ) . As noted by 
IFransson et al.l ( 20141 ). the NIR luminosity becomes domi¬ 
nant over the optical after the break in the light curve at 
- 370 d. 


3.3 Spectral evolution of SN 2010jl 

As seen in Fig. [T] the major features of the spectra of 
SN 2010jl remain mostly unchanged throughout its evolu¬ 
tion, with some notable exceptions. The continuum of each 
spectrum can be approximated by a T ~ 6000 — 8000 K 
blackbody spectrum, however, we note a slight bluing of the 
continuum at late times which may be due to contaminat¬ 
ing light from the host galaxy as the SN fades. We also 
detect [O ili] narrow emission features near 5000 A, which 
become more prominent at late times and are also likely 
due to contamination from the underlying star forming re¬ 
gion. Host galaxy spectra of two Hil regions ^ 7 arcsec 
north and south of the position of the SN from IStoh et ^ 
(1201 ih clearly show strong, narrow Balmer and [O III] emis¬ 
sion features indicative of recent star formation. A number 
of spectroscopic parameters, discussed below are given in 
Table HI 

The most prominent lines in the spectra are the strong, 
broad Balmer emission features. We detect the Ha through 
H(5 lines individually, while lines bluer than H(5 are blended 
together. Superimposed on the broad Balmer lines, we detect 
narrow blueshifted P Cygni absorption in the MODS spec¬ 
trum on t = 404 d and in the DIS spectrum on t = 512 d, 
but these components are not resolved in the rest of our 
lower resolution spectra. The width of the P Cygni compo¬ 
nent measured from the peak to the minimum is typically 
~ 100 km s“^. The evolution of the Ha prohle is discussed in 
detail in Section [3.3.II The broad He I AA5876 and 7065 lines 
also show evidence for a superimposed narrow component. 

The Balmer lines, particularly Ha, show a remarkable 
growth in strength relative to the continuum for t < 400 d 
as indicated by the equivalent width (EW) of the lines, mea¬ 
sured using the IRAF task splot. From t = 28 to 450 d, the 
magnitude of the EW of Ha strengthened from 167 zb 5 A to 
3400 zb 500 A , and it remained strong at ~ 2000 A in each 
of the later spectra (See Table H]). Similarly, the magnitude 
of the EW of H/3 strengthened from 35 zb 2 A to 300 zb 40 A 
in approximately the hrst 400 d. The increase in strength of 
the Balmer lines is indicative of interactions of the SN ejecta 
with a dense hydrogen rich CSM. We also note that the He I 
lines show similar evolution to that of the Balmer lines, but 
are weaker in strength. 

As given in Table 0] and shown in the left column of 
Eig.[5l we calculate the Ha luminosity, as well as the Ha/H^ 
and H/ 3 /H 7 luminosity ratios. Luminosities were calculated 
by subtracting a linear approximation to the continuum un¬ 
der the line, and then integrating the hux for the full width 
of the line. Eor t = 28,34, 512, and 775 d, the luminos¬ 
ity measurements are only approximate due to insufficient 
wavelength coverage in the spectra for absolute hux calibra- 
tion using avai l able p hotometry. We hnd similar results to 
lEransson et al.l (|2014l ) (compare with their hg. 12). The Ha 
luminosity increases by a factor of ~ 1.5 from t = 28 to 
154 d. After the break in the light curve, which occurred 
when the SN was close to the Sun, the luminosity drops 
again, and continues to fall by a factor of ~ 10 between 
t ~ 400 and t ~ 900 d. Erom the power law hts shown in the 
hgure and discussed above in Section [321 Ihe maximum in 
the Ha luminosity likely occurs at t ^ 350 d. This is ^ 10 d 
earlier than the date found by lEransson et al ] (l2ni4h . but 
the lack of data between t ~ 200 and t ~ 400 d makes it 
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Table 4. Spectroscopic parameters of SN 2010jl 


Phase 

EW (A) 

Hq H/3 

EW (A) 

Hel(5876 A) He I (7065 A) 

Ha luminosity 
(10^^ erg s“^) 

Hq/H/3 

H/ 3 /H 7 

28 

-167 (5) 

-35 (2) 

-14 (4) 


6.4 

4.0 

1.8 

34 


-40 (3) 





1.6 

55 

-310 (10) 

-59 (6) 

-24 (1) 

-11 (1) 

5.8 

5.2 

2.5 

81 

-480 5o) 

-80 (6) 

-41.6 (0.7) 

-18 (3) 

7.3 

5.2 

2.6 

154 

-950 5o) 

-134 (9) 

-54.6 0.8) 

-22 (2) 

9.6 

6.1 

3.1 

404 

-2700 5oO) 

-300 (40) 

-51 (2) 

-24 (1) 

7.3 

6.2 

3.6 

450 

-3400 5oO) 

-150 00) 

-75 (2) 

-36 (2) 

7.1 

14.3 

3.5 

468 

-2400 000) 

-270 00) 

-42 (1) 

-20 (1) 

5.6 

6.6 

4.0 

501 

-2000 000) 

-190 (30) 

-33.7 (0.9) 

-19 (2) 

4.3 

7.4 

3.8 

512 

-2400 000) 

-200 00) 



4.1 

8.0 

5.1 

775 

-1900 000) 

-98 (5) 

-18 (2) 


1.3 

10.8 

2.9 

824 

-2400 000) 

-126 (5) 

-24 (2) 

-18 (3) 

0.86 

12.9 

3.4 

891 

-2100 000) 

-52 (2) 

-11 (2) 

-12 (2) 

0.52 

26.1 

3.6 

909 

-3000 000) 

-150 (7) 

-34 (2) 

-14 (4) 

0.59 

16.3 

2.7 


The phase is the number of days since earliest detection on JD = 2,455,479.14. l-cr uncertainties are given in 
parentheses. The luminosities at t = 28, 34, 512 and 775 d are only approximate due to insufficient spectral 
coverage for flux calibration with synthetic photometry. 



t (days) t (days) 


Figure 5. The evolution of various features of the Balmer emis¬ 
sion lines. As in Fig. 3] open points are only approximate. Upper 
left: The integrated Ha luminosity and corresponding power law 
fits. Lower left: The Ha/H/3 and H/ 3 /H 7 luminosity ratios are 
shown as cyan squares and magenta dots, respectively. Upper 
right: The FWHM of the IW component of Ha. Lower right: The 
blueshift of the IW component of Ha measured as the median 
velocity of the IW profile. 


difficult to put strong constraints on this date. As shown 
in Fig. [5l the Ha/H/3 ratio increases from ~ 4.0 to ~ 16.3 
between t = 28 to 909 d. The H^/Hy ratio increases slightly 
from ~ 1.8 to ~ 5.1 between t = 28 and 512 d, and then de¬ 
creases slightly again to around 3.0 from t ~ 800 to ~ 900 d. 

An obvious change in the appearance of the spectra 
is the development of blended Call IR triplet emission 
(AA8498, 8452, 8662). These lines are absent in the earli¬ 
est spectra, and begin to develop by t = 81 d. The Call 
IR triplet is a more prominent feature in the spectra by 
t = 154 d. These emission features may be blended with a 


number of other lines including the O I A8446 feature. We 
also detect the presence of [Neiii] A3865.9 in each of our 
spectra with sufficient wavelength coverage. 


3.3.1 Evolution of the Ha emission profile 


The evolution of the Ha emission line is illustrated in Fig. [6] 
from t = 28 to t = 909 d. A summary of the spectroscopic 
parameters of the emission features, including Ha, is given 
in Table m As discussed in Section [3.31 a remarkable feature 
of the Ha emission is the increase in the EW of this line 
until t = 450 d, and the continued strength of the line with 
respect to the continuum throughout the rest of the spectra, 
indicating strong, successive interactions of the SN ejecta 
with a dense hydrogen-rich CSM. 

To examine the evolution of the Ha profile in detail, we 
performed multi-component fits to the profiles at each epoch 
using a Levenberg-Marquardt least-squares minimization 
algorithm. A linear approximation to the continuum near 
the line was subtracted from each profile before performing 
the fits. We find the profiles to be best fit by three sepa¬ 
rate components: an NW Gaussian, an IW Gaussian, and a 
broad or IW Lorentzian. The results of our fits are super¬ 
imposed on the data in Fig. [6l and the parameters for each 
component, including the full IW component (sum of the IW 
Gaussian and IW Lorentzian), are given in Table [5] Previ¬ 
ous work decomposed the Ha profile into a two-component 
profile consisting of an NW component with full width at 
half maximum (FWHM) of < 1000 km and an IW com- 
ponent with FWHM velocity ^ 2000 — 3500 km s~^ (e.g., 
Patat et ^ 2011 : ISmith et al.l l2011bl . I 2 OI 2 I : IZhang et al.l 
2 OI 2 I : Fransson et al. l2014l i . In this data set we find evidence 
for an NW and IW component, but we are unable to achieve 
satisfactory fits of the peaks and wings of the profiles simul¬ 
taneously using only two components. Instead, we find that 
the IW component is itself usually best fit by two separate 
components, an IW Gaussian and an IW Lorentzian. The 
need for the extra component is illustrated further by the 
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Table 5. Ho: profile component parameters 


Phase 

NW Gauss 

Ha component FWHM 
(km s“^) 

IW Gauss IW Lorentz 

IW Full 

Shift 
(km s“^) 

28 

90 

320 

2200 

2200 

-153 

55 

360 

1850 

6660 

2450 

-403 

81 

340 

1850 

6340 

2590 

-417 

154 

360 

1990 

5470 

3070 

-528 

404 

170 

3140 

2030 

2430 

-791 

450 

390 

2860 

1810 

2350 

-744 

468 

320 

2770 

1620 

2240 

-651 

501 

350 

2590 

1430 

2210 

-663 

512 

80 

2060 

1260 

2100 

-646 

775 

150 

2020 

970 

1750 

-531 

824 

360 

2020 

910 

1760 

-416 

891 

290 

2040 

890 

1590 

-470 

909 

280 

1910 

890 

1480 

-506 


The phase is the number of days since earliest detection on JD = 
2,455,479.14. The shift is measured as the median of the full IW com¬ 
ponent. 


emerging asymmetry of the IW component at late times (see 
Fig. [ 8 ] and discussion below). 

The center of the NW Gaussian component for each Ka 
profile is typically consistent with 0 zb 100 km s“^ relative 
to the rest frame of the host galaxy. In our highest reso¬ 
lution spectra, the LBT/MODS spectrum from t = 404 d 
and the APO/DIS spectrum from t = 512 d, we see that 
the NW component is actually a P Cygni component con¬ 
sisting both of emission and blueshifted absorption with a 
typical width of ~ 100 km s“^ measured from peak to min¬ 
imum. For these profiles, the NW Gaussian is unable to fit 
the absorption component and is only representative of the 
emission component of the narrow P Gygni profile. In the 
lower resolution spectra, the P Gygni profile is unresolved, 
and is fit by the NW Gaussian with a typical FWHM of 
~ 300 km s“^. At late times as the SN fades, the NW com¬ 
ponent may suffer from significant contamination from the 
underlying star forming region. 

In the earliest profile from t = 28 d, the NW component 
of the line is best fit by both Gaussian components, one with 
FWHM velocity of ~ 90 km s“^ and one with FWHM veloc¬ 
ity of ~ 320 km s“^. This may be due to a partially resolved, 
narrow P Gygni component, resulting in the observed bump 
on the blue side of the narrow peak. The IW component of 
the profile from t = 28 d can be approximated by a single 
IW Lorentzian with FWHM velocity of ~ 2200 km s~^. 

In the subsequent profiles, the IW component, com¬ 
posed of the IW Gaussian and IW Lorentzian, shows a 
number of important changes as the profile evolves. In 
the profiles from t = 55,81, and 154 d, the high ve¬ 
locity wings of the IW component are dominated by a 
broader (FWHM ~ 6000 km s“^) Lorentzian component. 
By t — 404 d the high velocity wings have vanished, and 
the IW Lorentzian component becomes considerably nar¬ 
rower (FWHM < 2000 km s“^). By t = 450 d, the full IW 
component begins to be more dominated by the Gaussian 
component (FWHM - 2000 - 3000 km s“^). The FWHM 
of the IW component is plotted in the upper right panel of 
Fig. El as a function of time. The IW component increases 


in FWHM to ~ 3070 km s“^ at t = 154 d, before begin¬ 
ning to decrease for t > 404 d toward ~ 1500 km s~^ at 
t = 909 d. We note the qualitative similarities in the evo¬ 
lution of the Ha IW component FWHM and the integrated 
Ha luminosity shown in the upper left panel of Fig. [S] 

The profile is mostly symmetric in the Hiltner/OSMOS 
spectrum on t = 28 d. By t ~ 50 d, however, the pro¬ 
files show a progressive enhancement toward the blue, while 
the red side appears suppressed. We also note a progressive 
blueshift in the profiles, measured as the median velocity of 
the full IW component. The blueshift of the IW component 
of Ha is shown in the bottom right panel of Fig. [5] The pro¬ 
file becomes significantly blueshifted with an IW component 
median velocity of ~ 400 km s~^ on t = 55 d, and becomes 
more blueshifted with a median velocity of ~ 800 km s“^ on 
t = 404 d. For t > 404 d, the IW component becomes less 
blueshifted again, with a median velocity of ~ 500 km s~^ 
in the last profile on t = 909 d. 

The change in shape of the Ha profile with time is il¬ 
lustrated in the top row of Fig. [7| In the top left panel, 
the profiles are scaled to match in the velocity range from 
— 1200 km to —1000 km s“^. By t = 55 d the profile 
shows signs of a relative deficit on the red side of the IW 
component, and this trend continues until t = 154 d. The 
profile from t — 404 d is shown in black and no longer ex¬ 
hibits high velocity wings on either the red or blue side. A 
similar trend is observed for the H/3 and H 7 emission line 
profiles, shown in the middle left and lower left panels, re¬ 
spectively. The late time evolution of the IW component Ha, 
shown in the top right panel of Fig. [T] changes significantly 
for t > 404 d. The profiles, now scaled to match near their 
peaks, continue to narrow with time, but exhibit a progres¬ 
sive suppression of flux toward the blue side of the peak. 
This results in profiles which appear asymmetric about the 
peak by t = 775 d. 

In Fig. [H we analyze the symmetric properties of the 
Ha profiles in detail. The continuum subtracted profiles are 
shown in black as in Fig. [G] but the NW Gaussian com¬ 
ponent has now been subtracted for clarity. The reflection 
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Figure 6. The continuum subtracted Ho; profiles are shown 
in black, where 0 km s“^ corresponds to the host galaxy rest 
frame wavelength of Ha. The components of each fit (a NW 
Gaussian, an IW Gaussian, and an IW Lorentzian) are shown 
in green, yellow, and blue, respectively. The best fit to the 
data is shown as the red curve. Each profile is shifted up by 
6 X 10“^"^ erg s“^ cm“^ from the one below, and the last 

four profiles have also been multiplied by the factor shown in 
the lower right corner along with phase for clarity. For the days 
indicated with asterisks, the fluxes are only approximate, corre¬ 
sponding to the open points in Fig. [4| 


of the profile about the peak of the IW component is then 
shown in orange for comparison so that for a symmetric 
profile the orange and black curves will line up with each 
other. The profiles display some asymmetry from t = 28 
to 154 d, where the the flux to the red side of the peak 
appears suppressed, but it is possible this is due to uncer¬ 
tainties in estimating and subtracting the continuum. This 
asymmetry continues in the day 404, 450, and 468 profiles, 
but appears less pronounced. The profile from t = 775 d, in 
contrast, shows a clear suppression of flux on the blue side 
of the peak. This asymmetry persists in the three remaining 
profiles, but again becomes noticeably less pronounced with 
time. 

A number of scenarios have been proposed to explain 
the evolution of the Ha profile, and we examine these pos¬ 
sible interpretations below in Section 


4 DISCUSSION AND PHYSICAL 
INTERPRETATION 

Here we discuss the physical implications of the observed 
properties of SN 2010jl presented above, and provide ad¬ 
ditional commentary on physical interpretations previously 
presented in the literature. 


Figure 7. Sequence of emission line profiles for Ha (upper), H/3 
(middle), and Hy (lower). Early time profiles are shown in the 
left column, and late time profiles are shown in the right column. 
Profiles are scaled to match in the velocity range indicated by the 
dashed vertical lines in each panel. The last four Hy profiles from 
t = 775 to 909 d were not included due to low S/N. 


4.1 The NW CSM lines 

We observe the presence of narrow spectral emission features 
in the Balmer series as well as in the optical He I AA5876 
and 7065 lines. We interpret this as emission from slow- 
moving gas in the CSM, photoionized by light from the SN. 
In the higher resolution spectra from t = 404 and 512 d, the 
narrow component of the Balmer lines exhibit blueshifted 
P Cygni absorption by the expanding foreground CSM gas. 
The width of the P Cygni component, measured from peak 
to minimum, is typically ~ 100 km s“^, which we take as an 
estimate of the CSM expansion velocity. T his is consistent 
with t he estimate of the CSM velocity by iFransson et al.l 
(|2014l) , and with the suggestion of an LBV star as a possible 
progenitor (|Smith et al.ll201li3 ). 


4.2 The IW lines 


As discussed above in Sections 13.31 and 13.3.11 the broad 
emission lines of H and He show a remarkable growth in 
strength until t ~ 450 d and clearly indicate strong, succes¬ 
sive interactions between the SN ejecta and dense, hydrogen- 
rich CSM. Similar observations of SN 20 1 Ojl have been dis¬ 
cussed at length i n the literature (e.g ., Smith et al. _ 20121 : 
IZhang et al.ll2012l : IFransson et al.]l20l4lGah et al. 20141 ). 

The IW component of the Balmer lines, especially Ha, 
have been of particular interest because they show some 
evidence for newly formed dust in the post-shock shell of 
the SN, but not ah lines of evidence are consistent with this 
interpretation. 
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Figure 8. The continuum subtracted Ha profiles after subtract¬ 
ing the NW component for clarity (black). The reflection of each 
profile about the IW component peak is shown in orange. Each 
profile is shifted up by 2 x 10“erg s“^ cm“^ from the 

one below, and the last four profiles have also been multiplied by 
the indicated factor for clarity. Days since the earliest detection 
are listed for each profile. For the days indicated with asterisks, 
fluxes are only approximate as in Fig. [6] Note that the profiles 
appear approximately symmetric until at least t = 512 d, after 
which they exhibit an enhancement in flux to the red. 


In our earliest spectrum from t = 28 d, the IW com¬ 
ponent of Ha is well approximated by a single Lorentzian 
with a FWHM velocity of ~ 2200 km s“^ (Fig- EJ), which 
is characteristic of the scattering of Ha photons by thermal 
electrons in an optically thick CS M. This resul t is con sistent 
with the early time findings of ISmith et al.l (|2012l ) . Soon 
after, however, the IW Balmer profiles begin to show sig¬ 
nificant attenuation of the redshifted flux and a progressive 
blueshift of the line center with time until at least t = 154 d. 
Based on similar findings. ISmith et al.l (|2012l ) proposed that 
dust formation in the cool, post-shock shell of the SN may 
be obscuring red- shifted pho t ons f rom interactions on the 
far side of the SN. lGall et al.l (|2014l ) observed similar atten¬ 
uation of redshifted photons and progressive blueshifts of 
the line centers in a number of IW lines for approximately 
the first 240 d. They also show that this blueshift is stronger 
at shorter wavelengths, w hich would be expec ted if dust for¬ 
mation were responsible. iMaeda et al.l (|2013l ) advocate for 
post-shock dust formation as the mechanism responsible for 
the profile shapes with similar findings at t ~ 500 d. This 
interpretation is dependent on the assumption, though, that 
early time profiles accurately represent the intrinsic, unat¬ 
tenuated line profiles throughout their evolution, which my 
not be the case. 

As is apparent in Fig. [71 the profile shape and evolu¬ 
tionary trend of the Balmer emission features has changed 


significantly by t = 404 d. The underlying broad component 
(FWHM velocity ~ 6000 km s“^) present in earlier spec¬ 
tra has vanished, and the profiles continue to narrow with 
time throughout the remaining spectra. Rather than pro¬ 
gressive attenuation of the red side of the profiles, we now 
see an increasing suppression of the flux to the blue of the 
IW component peak, resulting in profiles which appear very 
asymmetric about the peak velocity by t = 775 d (Fig. [8]). 
The profiles remain blueshifted (i.e., the majority of the line 
flux still falls to the blue of the host galaxy velocity) even 
though the degree of the blueshift appears to be decreas¬ 
ing for t < 400 d. These results are apparently inconsistent 
with the scenario of continued dust formation at late times, 
since one would expect a progressive blueshift as light from 
the receding side of the SN is increasingly blocked by newly 
for med dust. _ 

iFransson et al.l (|2014l ) claim that even as the Ha pro- 
hles become more blueshifted, they retain a symmetric shape 
characteristic of electron scattering. They invoke electron 
scattering as the origin of the prohle shape, with a macro¬ 
scopic bulk velocity toward the observer in the scattering 
medium to explain the blueshift of the IW component, and 
reject dust formation scenarios altogether. Although the line 
prohles appear at least approximately symmetric about the 
peak velocity for the hrst t ~ 500 d, the clear asymme¬ 
try present in the prohles after t = 775 d (see Fig. [8]) can¬ 
not be accounted for by electron scattering alone. The ap¬ 
pearance of asymmetric and shifted line prohles could be 
due to asymmetries in the exp losion or CSM geo metries. 
Spectropolarimetry results from IPatat et al.l (| 20 111 ) also in¬ 
dicated an asymmetric geometry is likely, but we do not 
explore this possibility further in this paper. 

The formation of emission line prohles is dependent on 
numerous factors, including multiple velocity components 
and emission/absorption mechanisms. Additionally, when 
examining scenarios involving obscuration by newly formed, 
or pre-existing dust, it is impossible to know the true line 
prohle for each epoch, and assumptions about intrinsic pro¬ 
hle shape must be made. For these reasons, although the line 
prohle evolution appears inconsistent with the post-shock 
dust formation scenario at late times, we hnd it difficult to 
rule out any scenario based solely on the shapes and evolu¬ 
tion of emission line prohles. Further evidence from the light 
curves of SN 2010jl regarding the presence of newly formed 
or pre-existing dust are discussed below in Section 14.31 


4.3 The optical and NIR light curves 

As is evident from Fig.[31 the optical light curve of SN 2010jl 
displays a clear two-stage evolution. After the peak in the 
optical luminosity at t ~ 20 d, the light curve shows a 
steady decline which is well ht by an Lopt(t) = 1.02 x 
10^^(t/100 d)“°'^^ erg s~^ power law decay for at least 
t < 200 d. At t ~ 400 d, the light curve shows a clear 
break and is now well ht by a much steeper power law decay 
given by Lopt(t) = 6.01 x 10^^(t/370 d)“^'®^ erg s“^. Based 
on the intersection of the power law hts, the break in the 
light curve occurs at t ~ 370 d. At essentially the same time, 
the near-IR luminosity rises rapidly so that the evolution of 
the bolometric luminosity changes little. A problem in the 
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models of SN 2010jl to date is that they do not provide a 
natural explanation of this transition. 

lOfek et ^ (|2014l l lacked information on the near-IR 
evolution and based their model of the bolometric light curve 
as a simple correction to the PTF R-band light curve. Thus, 
the break in the optical light curve was interpreted as a rapid 
drop in the luminosity being produced by the CSM interac¬ 
tions. They considered three possibilities. First, the shock 
speed, and hence the luminosity, will drop rapidly once the 
swept up CSM mass begins to exceed the ejected mass (the 
so-called ‘snow plow’ phase). Second, the shock could simply 
have reached the edge of the dense regions of the CSM. Fi¬ 
nally, the CSM density eventually drops to the point where 
cooling becomes inefficient (the slow cooling stage) and a 
decreasing fraction of the shock energy is converted to opti¬ 
cal photons. They rule out the last solution based on their 
estimate of the CSM densities, and prefer the first possibil¬ 
ity to the second. However, with the addition of the infor¬ 
mation on the near-IR emission, there is no dramatic drop 
in the bolometric luminosity to be explained, or one must 
introduce some independent process to explain the sudde n 
rise in the near-IR e mission see n by Andrews et al. I il 2 niih. 
iMagda et alJ ll2013l '). iGall et ^ (l2014l ^ and iFranssonetal] 
il2ni4i. 

iFransson et alJ (l20lA interpret the two evolutionary 
phenomena as being unrelated. The optical luminosity drops 
because the shock is reaching the e dge of the r egions with 
a dense CSM or, as proposed by lOfek et ^ (|2014l b the 
shock slows as the swept up CSM mass becomes compara¬ 
ble to the ejecta mass, while the near-IR emission is caused 
by an echo from pre-existing dust radiatively heated by 
light from the SN peak. The first problem with this in¬ 
terpretation is that it requires remarkable coincidences in 
timing and energetics because the near-IR emission has to 
rise just as the optical emissions drops while maintaining 
a smooth evolution of the b olometric luminosity. The sec¬ 
ond problem, as discussed bv IFransson et aP (|2014l b is that 
the timing requires the inner edge of the dust shell to be 
at Rin ~ 6 X 10^^ cm, which makes it very difficult to have 
dust grains hot enough to produce the near-IR emission. 
The observed temperatures are Td ^ 1500-2000 K while we 
estimate that a = O.l/xm grains would hav e temperatures 
of onl y Td = 9 50/800 K for graphite /silicate lfTaor Draind 
lll99.tl grains. IFransson et ahl (|2014| j argue that this can be 
solved by making the grains very small, but even for a grain 
size of a = 10“^/xm, we estimate dust temperatures of only 
Td = 1100/680 K for graphitic and silicate grains. Further¬ 
more, for grains created in a wind, or any ejecta, there is 
a strong correlation between the grain size and the allowed 
optical depth because both are directly depend ent on the 
densit y - in essence r oc a oc M. In the models of iKochanekI 
ll2014[) . a wind producing such small grains would then have 
a negligible optical depth at these large radii and so could 
not support the significant optical depth needed to produce 
the observed near-IR flux as an echo. Even if we are gener¬ 
ous, the original dust-forming wind would have had at most 
Tv — 1 and the optical depth of this large radius shell would 
be smaller by RformiRin — 0.01 to leave an optical depth 
of Tv — 0.01 or smaller that simply could not support such 
a luminous echo. 

The only natural way to explain a sharp break in the 
optical light curve combined with a simultaneous rise of the 


near-IR emission to an almost identical luminosity is dust 
forming and becoming optically thick. It not only explains 
the simultaneity and the similar luminosities, but also the 
high dust temperatures because the temperature of newly 
forming dust is always clo se to the evapora tion temperature 
of 1500-2000 K. Following iKochanekI (|2014| j we can consider 
dust formation simply in terms of the ‘blackbody’ dust tem¬ 
perature aT^ij = L/IGttt^ set by the radiative flux. We will 
scale the results assuming dust forms at T^b = lOOOTbbs K 
and discuss Tbb and the actual dust temperatures in more 
detail below. We will model the bolometric luminosity evo¬ 
lution as L = 10 ^^L 43 t^QQ^^ erg s“^ where t = SOOtaoo d is 
scaled to the approximate time of the break in the optical 
luminosity evolution. 

As noted in other studies (e.g.. lAndrews et al.l l201ll : 
I Gall et ^ I 2 OI 4 I : IFransson et al.1 l2014h , dust cannot form in 
the material ejected during the transient because the dust 
temperature would be too high to form inside a radius 

i? c. 6 X cm (1) 

that is much larger than the ejecta radius 

R ^ 5 X 10^^x;e2t3oo cm (2) 


until long after 300 days if it is expanding at Ve = 
2000x;e2 km s~^. The evolution of these two radii are shown 
in Fig. [9] and the ejecta would not reach a large enough 
radius to allow dust formation until t ~ 2000 d. 

As discussed in iKochane'^ (|201l| j , the re-obscuration of 
the transients SN 2008S and NGC 300 OT-2008 can be ex¬ 
plained by dust re-condensing in their pre-existing dense 
winds. The GSM densities invoked to explain SN 2010jl are 
far denser than (even) the extreme AGB star winds invoked 
for these transients, although this is balanced by the far 
higher luminosities. Still, the optical and near-IR evolution 
of SN 2010jl strongly suggests that a similar process has 
occurred in SN 2010jl but under more extreme conditions. 
Inverting Equation [T] the time at which dust can begin to 
form at a given radius is 

tTo^ = O.ULlsR^fTd,! (3) 

where R = lO^^Rir cm. This is a very rapidly faffing func¬ 
tion of radius because the luminosity is dropping so slowly. 
In a wind with local density pw where we can ignore the 
outward expansion on the growth time scale ^ a grain 
can collisionally grow to radius a = fvcpwt^^'^^/^Pbuik 
where / = O.OI /2 is the condensible mass fraction, Vc — 
Ivci km s“^ is the collision velocity and pb ^ 2.2 g cm“^ is 
the bulk density of the grains. If we define pw = M/^.'KVyjR? 
by a pre-existing wind characterized by a mass loss rate 
M = MiM© yr“^ and wind velocity Vw = 100vw2 km s“^, 
a grain grows to radius 


a = 0.015 


Vw2Riy 


pm. 


(4) 


Very small grains have negligible visual opacities, so we can 
assume there is no significant absorption until a > amin — 
O.OOlttms A^m. This then gives the required growth time at 
any given radius. 


.grow 

^300 


0 gg 

f2MiVcl 


(5) 


which is a rapidly rising function of radius. On small scales 
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Figure 9. Scales related to dust formation. The steep, descend¬ 
ing dashed line shows the time dependence of the radius at which 
the radiation will permit dust formation assuming T 553 = 1 and 
setting L 43 = 0.64 to match our estimate at ^300 = 1 from Sec¬ 
tion [321 The ejecta radius (rising dotted line) does not enter the 
region where dust formation is possible until long after the ob¬ 
served time of the optical break at 370 d (horizontal solid line). 
For th e mass loss rate of Mi = 0.27vw2 invoked by lOfek et ah! 
il2ni4l'i . dust grains cannot have grown large enough to have a 
significant visual opacity below the second dashed line. At any 
given radius, dust will have a significant visual opacity only at 
times later than the sum of these two limits, which is given by 
the heavy solid curve. The resulting minimum time for dust for¬ 
mation is very close to the observed time scale of the break in the 
optical light curve and the dramatic rise in the near-IR emission. 


the radiation field keeps dust from forming and on large 
scales it simply grows too slowly. Fig.[9]shows these two time 
scales, t'^QQ and , after setting L 43 = 0.64 to match our 
estimate from Sec tion and Ml = 0.27vw2 to match the 
wind density from lOfek et al.l (|2014l i. 

Thus, there is a significant dust opacity at any given 
radius only once t > and 

Fig. [ 9 ] also shows the sum of these two time scales. Because 
of the different radial scalings of and there is a 

minimum time for dust formation. Moreover, without any 
particular manipulation of the parameters, this minimum 
time almost exactly matches the observed time of the optical 
break. If we minimize t/orm, we find that the radius at which 
dust first reaches size amin is 


R 


min 

17 


0.85 




^cl 


l/6^4/3 1/6 
^m3-^bb3 ^w2 


(6) 


and this occurs at time 


,min 

^300 


0.73- 


2/3 r 2/3 2/3 
(^rn3W3 ^w2 

2/3,V2/3^8/3 2/3' 


r^M- 


'■663 '^cl 


(7) 


If we identify this time scale with the break in the light 
curve, tb = then the required mass loss rate of 


Figure 10. As in Fig. [9] but with the total luminosity breaking 
from L oc to oc Shortly after the luminosity break, 

dust formation will occur over a broad range of radii. 


M - 0.62 y^-1 

hvciVMil 

is quite close to the values of M ^ 0 . 27vw2 and 

^ 0.11vw2 yr~^ invoked by lOfek et ^ (|2014l ) and 
iFransson et al ] (l2014h to explain the observed luminosity. 
Once there are grains of significant size a > amin, the op¬ 
tical depth becomes significant very quickly and there will 
always be hot, newly forming dust to support the near-IR 
emission. 

We should note that if a den se C SM exists on these ra - 
dial scales, the lOfek et al.l (|2014l ) and IFransson et ^ (|2014l ) 
scenarios for the origin of the break make dust formation 
near the time of the break even more inevitable. If either 
the total luminosity or the radiation temperature (or both) 
abruptly decrease, the radius outside of which dust forma¬ 
tion is permitted drops precipitously compared to the sce¬ 
nario shown in Fig. (9] At these smaller radii the CSM densi¬ 
ties are higher and the grain growth times are so short that 
dust formation occurs almost instantaneously once permit¬ 
ted by the radiation field. For example. Fig. [TO] shows the 
effect of holding Tbb fixed but adding a break in the lumi¬ 
nosity from roughly L oc to L oc we found in Sec¬ 

tion (321 Shortly after the luminosity break, dust formation 
occurs over a broad range of radii by the time t ^ 1000 d. 

While our solution naturally explains the observations, 
its parameters are somewhat strained. First, given the high 
radiation temperature (8000-9000 K) a dust formation tem¬ 
perature of Tbb — 500 K seems more likely than the Tbb = 
1000 K we have used in the scaling solutions (see iKochan^ 
I 2 OI 4 I ) . But lowering Tbb by a factor of two would require rais¬ 
ing M by a factor of 16 to keep the fixed. Note that the 
dust emission would still be very hot even at Tbb — 500 K 
because the smallest forming grains would be stochastically 
heated to temperatures close to their evaporation tempera- 
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tures and would radi ate most of their energy in these tem¬ 
perature spikes fsee iKochanekI [2014l h Second, once dust 
forms, the optical depth becomes very high, with 


Tv — 600 




/•1/6 1/6 5/6 ^1/3 
12 ^cl 


3 7-1 

^43 


(9) 


where we have scaled the visual opacity to nv — 
100 /^V 2 cm^ g“^ and this would be the optical depth for 
a wind extending from outwards. If it forms in a fi¬ 

nite annulus, the optical depth is reduced by 1 — RiniRout 
which will quickly become a modest correction (see Fig. [9]). 
Empirically, the observed optical depths are far smaller. 

A plausible solution is that the CSM has significant 
density inhomogeneities. For example, if we drop the re¬ 
quired Thb to 500 K and keep M fixed while putting fraction 
F = 0.01 of the CSM in clumps with an over-density of 
^ = 16, then dust will only form in the dense clumps and 
the mean optical depth would drop by the factor F. The 
effective drop would be still larger because of photons escap¬ 
ing by scattering though the low density regions. Combined 
with the possibilities allowed by the global asymmetries in¬ 
voked by other studies of this system, this provides a means 
of making the properties of our solution less extreme. 


5 SUMMARY AND CONCLUSIONS 

The luminous and long-lived Type Iln SN 2010jl shows 
strong evidence for the interaction of the SN ejecta with 
a dense, hydrogen-rich CSM shell. We obtained 12 optical 
spectra between t = 55 and 909 d, which adds significant 
coverage to the spectroscopic time sequence for this event. 
We also obtained an epoch of late time BVRI photometry, 
and derived a semi-bolometric, optical light curve through 
t = 775 d. 

Analysis of the spectra revealed NW (FWHM ~ 
100 km s“^) emission features of H and He associated with 
the slowing expanding CSM shell, ionized by the initial flash 
of the SN. The IW (FWHM ~ 1000 — 2000 km s“^) compo¬ 
nent of the Balmer lines also showed complex evolution, par¬ 
ticularly a progressive blueshift for at least the first ~ 150 d, 
and a pronounced asymmetry after ~ 750 d. These observa¬ 
tions are not fully consistent with any emission and broaden¬ 
ing mechanisms yet proposed, including post-shock dust for¬ 
mation and electron scattering in an optically thick medium. 

The optical light curve exhibits a steep drop off for 
t > 300 d that has been interpreted as a transition of the SN 
to the snow-plow phase, or as the emergence of the shock 
from the CSM. However, in light of the timing coincidence of 
the break in the optical and rise in the NIR light curves, the 
most natural interpretation is that the optical light from 
continued CSM interactions begins to be reprocessed into 
the NIR after the break. Using a simple collisional model 
for dust growth, we find that dust naturally reforms in the 
dense, pre-existing CSM wind on the observed timescale of 
the breaks in the optical and near-IR light curves with¬ 
out any fine tuning of the parameters. The model would 
be improved if the CSM is significantly clumpy, consistent 
with suggestions that it has more compl e x structure (e.g. , 
Chandr^j^jL 1 201 2d : iMaeda et al.l 1 20131 : lOfek et ^ I20l4 
Fransson et al. 2014 . We also emphasize that this inter¬ 


pretation indicates continued strong interactions of the SN 
ejecta and the CSM for > 2.5 yr. 

It is also interesting that a number of observational 
changes in the evolution of SN 2010jl are approximately co¬ 
incident with the optical break and rise in the NIR light, 
including the decline in the Ha luminosity, the loss of a 
broad, high-velocity component in the Balmer emission pro¬ 
files and subsequent narrowing of the lines, and the recession 
in the blueshift of Ha. Despite the wealth of studies in the 
literature dedicated to this fascinating object, the details 
of the true nature of SN 2010jl’s evolution remains contro¬ 
versial, emphasizing the importance of multiwavelength and 
continuous monitoring of such events. 
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